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Abstract 

We consider a quantum graph where the operator contains a po- 
tential. We show that this operator admits a heat kernel. Under some 
assumptions on the potential, this heat kernel admits an asymptotic 
expansion at t = with coefficients that depend on the potential in a 
universal way. These coefficients are spectral invariants, we compute 
the first few of them. 

1 Introduction 

A metric graph is a combinatorial graph where each edge is equipped with a 
length. A pair of a metric graph together with a differential operator acting 
on it, is called a quantum graph. The operator is of the form Laplacian plus 
lower order parts on each individual edge. Some boundary conditions are 
imposed at the vertices to make the operator self adjoint. 

Quantum graphs were introduced in the 1930s in the physics literature, 
they are a popular model for various processes involving wave propagation. 
If the operator has no potential, the eigenvalue equation on each edge can 
be solved explicitly, the eigenf unctions are sine waves. In this setting one 
can write down an exact trace formula. It relates the set of eigenvalues 
to topological properties of the graph. The first such trace formula was 
found by |Rot84] . since then various generalizations have been shown, see 
|KS99bl IKPS07I IBE09] or |BE08j for a survey. One can get an exact trace 
formula for an operator with potential, however it uses the abstract existence 
of solutions to the eigenvalue equation on each edge as input |RS12j . Thus 
without specific knowledge of what these solutions look like, one cannot read 
off any information about the graph or the potential from this trace formula. 

For this reason we consider the heat kernel. Let G be a quantum graph 

with the operator D := — + U{x)j acting on it, here f/ is a real potential 
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function. Non-Robin boundary conditions are imposed at all vertices, see 
definition 12.11 The heat kernel e{t, x, y) is a function that satisfies 



where 6 is the Dirac-5-distribution. We will study its asymptotic expansion 
for time close to zero. 

In a setting without potential, this approach has the immediate disadvant- 
age of being only asymptotic compared to the exact trace formula. However, 
as we will show, the coefficients in the expansion depend on the metric graph 
and the potential in a universal way and can be computed explicitly. This 
means they give rise to a series of spectral invariants such as the total edge 
length of the graph, the integral of the potential over the graph, as well as 
the sum of the values at the vertices, weighted with some factor depending 
on the boundary conditions. 

More specifically we get the following two theorems. 

Theorem 1.1. 

Assume the potential function U is in L^(G', M), then the operator D admits 
a heat kernel e{t, x, y) and this heat kernel is unique. 

Theorem 1.2. 

Assume that U is smooth on the edges, that U^'^''^^\v) = for all vertices 
V G V and I E N and that U'^'^^\x) is a continuous function on the entire 
graph G for all I G N, then this heat kernel admits an asymptotic expansion 
of the form 



Here the a^°' come from the boundary conditions, see \2.1\ The coefficients 
an{x) and a'^{v) are local and universal. With our assumptions on the po- 
tential, a^(f) = whenever n is not an integer. The first few non-zero 
coefficients are as follows. 






.n 
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This situation closely mirrors the case of Riemannian manifolds where 
similar questions have been studied for quite some time. One of the earliest 
results is the work by Minakshisundaram and Pleijel in the 1950s, |MP49t 
IMin53] . They show the existence and universality of the asymptotic expan- 
sion for Riemannian manifolds without boundary and compute the coeffi- 
cients up to 02 for the standard Laplacian, see also |MS67t IBGM71] . Their 
approach uses a parametrix, this is the easiest way to start out but has ser- 
ious drawbacks if one wants to consider manifolds with boundary and the 
computation of more coefficients gets prohibitively complex. 

Existence of a heat kernel and an asymptotic expansion with universal 
coefficients has been shown for Riemannian manifolds with boundary, arbit- 
rary elliptic operators and an extensive set of possible boundary conditions 
in |Gre71j using a method based on the inversion of the operator on the level 
of symbols. Using this broad existence result, |Gil79] computed the coeffi- 
cients up to as and ag for manifolds with boundary and operators involving 
a potential using recursion relations. 

We will first show existence of the heat kernel for a quantum graph using 
properties of the eigenf unctions. While this method is very general, it seems 
unsuitable to study further properties of the heat kernel. The heat kernel is 
a local object, something that is rather difficult to see from the eigenfunction 
construction as eigenf unctions and eigenvalues are highly non-local. 

We will therefore use a different technique, namely a parametrix. This 
comes at the cost of some restrictions on the allowable potential functions but 
in return gives us much more detailed information about the heat kernel. Us- 
ing a parametrix makes it obvious that the heat kernel is a local object. The 
parametrix approach also gives us the existence of the asymptotic expansion 
essentially for free and provides a method to compute the coefficients. 

The construction of a parametrix for the heat kernel on a quantum graph 
is done in three steps. First, we look at the construction on the real line. 
Second, we use the construction on the real line to build a parametrix on 
a star graph, ie multiple half-infinite edges with a common central vertex. 
Finally, we use a partition of unity argument to construct the parametrix for 
a graph. 

Once we have build a parametrix for the graph we can build a heat kernel 
from it with a method essentially analogous to the manifold setting. The 
parametrix approximates the heat kernel, so we can use it to show the exist- 
ence of an asymptotic expansion and compute the first few coefficients. 

This paper is organized as follows. Section [2] contains the general setup 
and the definitions of the boundary conditions we are going to use. Section 
[3] contains the general existence proof of the heat kernel. Next, section H] 
contains the construction of a parametrix in the three steps outlined above. 
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In section [5] we build the heat kernel for the graph from a parametrix. Finally, 
in section |6] we show the existence of an asymptotic expansion and compute 
the coefficients. 



2 Setup 

Let G be a metric graph. The number of edges and vertices is finite and 

the length of each edge is finite. Let the operator D := — + U{x)^ 

act on it, here [/ is a potential function. At this stage we only assume 
U G L^(G,R) and make no assumptions on its behaviour at the vertices. 
The sign convention for U might seem unusual but it has the advantage of 
avoiding any minus-signs in the coefficients in the heat asymptotics. 

We will set up boundary conditions at the vertices as in |KS99aj and make 
use of their classification. 

Definition 2.1. 

Given a vertex v of the quantum graph and a function f on the graph, 
enumerate the edges adjacent to v from l,...,a. Denote by fy the vector 
{fi{v), . . . , fa{v)) and by fj the vector {f[{v), . . . , fa{v)) where all derivat- 
ives are oriented away from the vertex. Then the boundary conditions can be 
written as 

■^v fv + /?; = 

for two matrices and B^. Let 

a"" := -(A + ikB^)-\A^ - ikB^) 

for A; G M. // the matrix a"" is k-independent, the boundary conditions are 
said to be non-Robin. This means, they don't mix conditions on the function 
with conditions on its derivative. Note that (t"" is unitary and that (a^)^ = Id, 
FKWOT^ . 



Collecting the matrices Ay and B^ for all the vertices as blocks along the 
diagonal in the matrices A and B, the operator D is self-adjoint if and only 
if AB* is self adjoint and {A, B) has full rank. 

We will henceforth assume that non-Robin boundary conditions are im- 
posed at all vertices. 

Example 2.2. The Kirchhoff- Neumann boundary conditions at a vertex cor- 
respond to the matrix 

2 

deg(^;) 

where 6aj3 is the Kronecker-(^. 
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Remark 2.3. One could consider a more general operator of the form 

Da-= (^i^ + A{x)^ -U{x) 

with an additional magnetic potential A but this magnetic potential can be 
completely absorbed in the boundary conditions through a gauge transform- 
ation, see |KS03] . 

We have = U^^DqU where U is unitary. The boundary conditions are 
transformed as A AU and B BU^ in particular if Dj, has non-Robin 
boundary conditions at all vertices then so has Dq. 

The eigenfunctions /_4 of the operator can be obtained from the ones 
of Dq by multiplying with a prefactor 

Ux) = e-'^^^^^'^-'^'^'Ux) 

Thus the only change in the asymptotics of the heat kernel would be this 
prefactor, we will therefore only consider operators without magnetic poten- 
tial. 

Definition 2.4. 

We write f{t) = 0{t^) if there exist constants T and C such that 

\f{t)\^ct^ vtG [o,r] 

We will also write f{t) = if f{t) = 0{t^) for all k ^ 0. Note that 

t^ = 0{t^) if and only if k' ^ k. If the function f depends on multiple 
variables we assume C and T are independent of the other variables. 



3 Existence of the heat kernel 

In this section we will prove existence of the heat kernel in a very general 
setting, the proof technique is specific to quantum graphs, it works because 
so much is known about the eigenvalue equation on the unit interval. 

Definition 3.1. 

Let {k'^,dj}j be a spectral resolution of D. The kj are the eigenvalues and 
the eigenfunctions 6j form an orthonormal basis o/L^(G,R). 

The existence of such a spectral resolution is well known, it follows from 
the fact that D is elliptic and the graph is finite, the proof is written out 
in |Kuc04] for the case of an operator without potential but it immediately 
generalizes to our setting. 
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Consider an individual edge without boundary conditions, parametrized 
as the interval [0, L], then the equation 



y"{x) - U{x)y{x) = P 



y{x) 



has two real linearly independent solutions and 2/2 fo^^ ^iny fixed value of k, 
|PT87] ■ If the potential is only in L^[0, L] these solutions are in C^[0, L] with 
an absolutely continuous derivative but not necessarily in C^[0, L]. Normalize 
them so that 



The eigenfunctions on each edge of the graph are then linear combinations 
of these two solutions, which we write as 



for some real parameters pj and aj. To get an eigenfunction on the entire 
graph one needs to pick the coefficients pj and aj on each edge in such a way 
that the boundary conditions are satisfied at all vertices. 

Remark 3.3. The construction above reduces the problem of finding eigen- 
values and eigenfunctions to a finite list of linear equations. This, or a very 
similar approach, is the basis of the proof of an exact trace formula for 
quantum graphs, |KS99bl IBMM [RST2] . 

Theorem 3.4. 

There exists a constant C independent of x (z G and j such that 



In other words, there exists a global bound on the size of all eigenfunctions 
of a quantum graph. 




Theorem 3.2. 



PT8V The two solutions y\ and can be written as 




^ji^) = Pj [cos{aj)y^' (x) + sm{aj)y2' (x) 



e,{x)\<C 
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Proof. As we assumed the eigenfunctions to be orthonormal we have 



\9j{x)\'^dx ^ 



G 



This imphes 



p1 



cos{aj)y'^^ (x) + sm{aj)y2^ (x) ) dx ^ 1 



f (cos{kjX — aj) + cos(aj)-Ef^ (x) + sm{aj)E2' (x)) dx ^ — 
Jo ^ ^ Pj 

11 

/ cos^(A;,x — a,) + 0(—)dx ^ 
Jo ' ' k 



In particular for k sufficiently large 



Pj 



As the kj are discrete, this implies there is a global bound for the p^'s. By 
theorem 13.21 there exists a constant C", independent of x and k, such that 

\y'i{x)\<c' \y'2{^)\<c' 

This implies the claimed bound. □ 
Theorem 3.5 |Gri07] Weyl law. 

Let G be a quantum graph with non-Robin type boundary conditions at all ver- 
tices and spectrum {kj}j. Then one can estimate the number of eigenvalues 
in any interval {Kq, Ki) by 



i^{j\Ko<k,<Ki}--{Ki-Ko) 

71 



< 2E 



where L is the total edge length of the graph and E is the number of edges. 



Theorem 3.6. 

Let G be a metric graph with the operator D :- 



+ U{x)] acting on 



it. Let {kj,6j}j be a spectral resolution of D. Then the heat kernel for D 
G exists, is unique and is given by 

oo 

e{t,x,y):=J2^~'^%i^WM 

j=0 



on 
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Note that if the potential is only assumed to he in L^(G', M) then the heat 
kernel is not necessarily smooth in the x and y variables. 



Proof. This is a well known formal expression for the heat kernel. The pre- 
vious two theorems imply that the infinite sum converges absolutely for all 
t > and aX\ x,y E G. As the eigenf unctions are smooth on the edges and 
satisfy the boundary conditions at the vertices so does the heat kernel. 

To get uniqueness we argue as in |BGM71j . If e{t,x,y) is a heat kernel, 
it can be written as 

oo „ 

eit,x,y) = y2fj{x,t)ej{y) with fj{x,t) = e{t,x,y)9j{y)dy 
These coefficients fj satisfy 



G 



dtfjix,t)=l dte{t,x,y)ej{y)dy 

Dye{t,x,y)ej{y)dy 
- I e{t,x,y)Dyej{y)dy 



G 



G 

= — kjfj(x, t) 

where we used the fact that D is self-adjoint. This implies 

fj{x,t) = g,{x)e-''^' 

Considering the limit t — )■ now shows gj{x) = Oj{x) and establishes unique- 
ness. □ 



4 Construction of a parametrix 

We will construct a parametrix in three steps, first on the real line, then on 
a star graph and finally on the quantum graph. 

Definition 4.1. 

A parametrix is a function hk{t, x, y) that satisfies the following: 

1. It is smooth for t > . 

2. It is smooth for t ^ if x ^ y and satisfies hk{t,x,y) = 0{t°°) for y 
bounded away from x. 
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3. {dt+Dy)hk{t, x,y) is smooth fort ^ and satisfies {dt+Dy)hk{t,x,y) = 
0{t^~'^). 

4. limt^o hk{t, X, y) = 5x{y) 

We will explain in the respective sections what we mean by smoothness 
on the star graph and the quantum graph. 

4.1 The real line 

We will start the construction of a heat parametrix on the real line. Let 
f{t,x,y) := -yi^e 4^ denote the heat kernel of the standard Laplacian 
-dl on R. 

We will use the ansatz 

k 

hkit, X, y) := fit, X, y) ^ Ui{x, y)t^ (1) 

1=0 

where the -u-functions are determined recursively. We will always assume 

A; ^ 1, SO there is at least one term beyond :^7j^e 4t . 
We have 



{dt-dl-U{y))hk{t,x,y) 



y 

(k-l k 
«=0 «=0 1=0 1=0 

We order the terms by powers of t and set all but the highest power to 
zero: 

={x - y)dyUo 

=(/ + + {y - x)dyUi+i - dlui - U{y)ui 

— dyUk — U{y)uk is the highest order term 

This leads to UQ{x,y) = const, we will show below in lemma that 
UQ{x,y) = 1. For / ^ 1 this leads to the transport equation. 

Lemma 4.2. 

The transport equation on the real line is 

= lui{x, y) + {y- x)dyUi{x, y) + DyUi^i{x, y) 
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where I ^ 1 and Uo{x,y) = 1. It has the solution 



y 



ui{x,y) = -{y - x) ' / {z-xY D^ui^i{x, z)dz 



If the potential U is smooth then so are all the u-functions. 
Proof. The homogeneous equation 

= lui{x, y) + {y - x)dyUi{x, y) 
can be solved with separation of variables. It has the solution 

u1^^{x,y) = c{y-x)-' 

The inhomogeneous equation is then solved using variation of constants. The 
smoothness follows from lemma 14. 5[ below. □ 

Remark 4.3. The general inhomogeneous solution consists of the particular 
one above plus a solution to the homogeneous equation. As we want our 
solution to be well defined for x = y we have to choose the constant c = in 
the homogeneous equation. Thus the given solution is the only solution for 
the ui that is relevant to us. 

Remark 4.4. The solution to the transport equation implies that Ui{x, y) only 
depends on the potential in some small neighbourhood around the interval 
from X to y, in particular ui{x, x) only depends on a small neighbourhood of 
X. We will refer to this property by saying that the w-functions are local. 

Lemma 4.5. 

Let f be a function, and let 

g{y) = {y-a)-^ !\z - af-' f{z)dz 



then g is at a and 

9{o) 



f{o) 



n 

Proof. Use the e — 6 criterion for the C° case and a Taylor expansion for 
higher orders. □ 

Lemma 4.6. 

Let hk{t,x,y) be defined as in equation ^ with the u-functions defined by 
the transport equation, lemma \4^ Assume the potential U is smooth. Then 
hk{t,x,y) is a parametrix on the real line. 
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Proof. The first three properties just follow from the construction and the 
smoothness of the ui. 

For the last one recall that /(t, x, y) = -^^^e 4t is the heat kernel for 
the Laplacian without potential. Let A^^^ be an open neighbourhood of x with 
compact closure. Then 




= 'iIj{x)uq{x, x) 



as each integral in the sum converges to ui{x,x)'i/j{x) by properties of /, so 
we have to set uo{x, x) = 1 . □ 

4.2 A star graph 

A star graph consists of d half-infinite edges glued together at a central vertex 
V. At this vertex non-Robin boundary conditions are imposed as in section 

m 

Definition 4.7. 

We say a function is smooth on the star graph if its restriction to any two 
edges gives rise to a smooth function. 

This implies that all odd derivatives of the function vanish at the central 
vertex, whereas the even derivatives are continuous at the vertex. 

Assumption 4.8. 

We will assume that the potential function U is smooth on the star graph. 

Remark 4.9. Through carefully keeping track of degrees of differentiability 
the assumptions on the potential can be relaxed. In this case, only the lower 
order w-functions are well defined and they have low degrees of differentiab- 
ility at zero. However, any approach that uses a parametrix will require at 
the very least that the potential is continuous and that its first derivative 
vanishes at all the vertices. The general proof of existence of the heat kernel 
in section |3] indicates that even these conditions should not be necessary. 
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We will therefore not try to find the weakest possible differentiability 
conditions for the parametrix approach but keep the assumption above. This 
will significantly simplify further computations as this assumption makes all 
the w-functions smooth. Hopefully some different approach to heat kernel 
asymptotics on quantum graphs will give rise to results without such strong 
conditions on the potential in the future. 

Let Xa denote the coordinate for the edge a, parametrized as [0, oo). 

In order to model a reflection at the vertex we will continue the potential 
function via Ua{—Xa) '■= Ua{xa) onto the entire real line. The assumptions 
on the potential make this a smooth continuation. 

Let f\^(t, Xa, —yp) denote the parametrix for the real line for the potential 
Ua{x) for X ^ and Up{—x) for x ^ 0. 

We define the parametrix of a star graph as follows. 

Definition 4.10. 



where b^j^ is the Kronecker-6 and a'^^ is an element of the matrix with the 
boundary conditions, see definition \2.1[ This implies that satisfies the 
boundary conditions at the central vertex. 

Lemma 4.11. 

The function hl.{t,x,y) is a parametrix on the star graph. 

Proof. The first three properties follow from lemma 14.61 For the last one, 
assume that x is not the central vertex and lies on the edge a. Let Ga denote 
the restriction of the graph to the half edge a and A^^ an open neighbourhood 
of the point x with compact closure. Then 



The first term converges to '^(xq,) and the remaining ones all converge to 
zero by lemma 14.61 The case of x equal to the central vertex follows by 



Let 



hl{t,Xo,,yi3) : 



^a/sK^it, Xa, ya) + (T^^hf{t, X^, - J//?) 




continuity. 



□ 
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4.3 Construction of a parametrix for a quantum graph 



We assume that the graph does not have loops or multiple edges. This can be 
achieved without loss of generality by inserting additional vertices of degree 2 
with Kirchhoff-Neumann boundary conditions. These will not influence the 
spectrum of the operator so they will not show up in the asymptotics but 
they make the definition of the parametrix easier. 

Let Iq denote the length of the shortest edge on the graph and let d denote 
the distance function on the graph. 

Definition 4.12. 

Let ?7 : M — )■ [0, 1] be a smooth cut-off function such that 

• V l(-oo,i]= 1 l[|,oo)=0 

• — x) = 1 — rj{x), that is rj is symmetric about the point (|, |) 
Definition 4.13. 

Let := {x E G \ d{x,t)) ^ ^} be a neighbourhood of the vertex v. Let 

3 

Vio '■= U^gyV;^ and denote the complement of this set in G by Vi^. 
33 3 

Definition 4.14. 

We define a partition of unity {Xv}vev on the graph G as follows. For x on 
an edge adjacent to v set 

Xv{x) := r]{l{x)~^d{x,v)) 

where l{x) is the length of that edge. Set XvW) — ^^v' f^f^ vertices and 
Xv{x) = if X is not on an edge adjacent to v. This implies Xv Iv = 1- 

For the potential function, we will just carry over the smoothness assump- 
tion from the star graph to an arbitrary graph. 

Assumption 4.15. 

We assume that the potential function U is smoooth on the graph. In partic- 
ular 

[/f+i)(^;) = and U^^'\v) = U^J\v) 
for all edges a, a' adjacent to v for all vertices v, for all I e N. 

This guarantees that the w-functions are well defined and their restrictions 
to any two adjacent edges are smooth functions. 
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Definition 4.16. 

Let x,y E G, then we set 



hkit, X, y) := 77(2/0 ^d{x, y)) ^ Xv{.x)hl{t, x, y) 



where hi is the parametrix on the star graph constructed in definition 4J_0. 

Remark 4.17. Note that this is well defined in the sense that the factor in 
front of hl{t, X, y) is only non-zero if both x and y lie on edges adjacent to v. 

The notation here is global and coordinate free. The coordinates only 
come into play at the level of the star graphs. If a point lies on an edge 
between the vertices v and v' this edge will be parametrised from v to v' 
in hi but from v' to v in h^. so the same point has different coordinates in 
different elements of the sum. 

Lemma 4.18. 

The function satisfies the boundary conditions at all the vertices. 

Proof. If d{x,y) > y then hk{t,x,y) = because of the initial cut-off func- 
tion. If d{x, ^ ^ and x and y are both within a V^^ neighbourhood 

"3" 

of a vertex v we have hk{t,x,y) = hl{t,x,y) which satisfies the boundary 
conditions by construction, see definition 14.101 □ 

Lemma 4.19. 

The function hk{t,x,y) is a parametrix on the graph G. 

Proof. The first two properties follow directly from the analogous properties 
of hi in lemma 14.111 For the third one we have 

{dt + Dy)hk(t,x,y) 
=ri{2lQ^d{x, y)) ^ Xvix){dt + Dy)hl{t, x, y) 



dy{r]{2lo y))) ^ Xvix)dyhl{t, x, y) 
9j(r/(2/o-M(x, y))) ^ Xv{x)hl{t, x, y) 



The first term extends smoothly to t — > and satisfies {dt + Dy)hl{t, x, y) = 
0{t^~'i) by lemma H. Ill The second and third term are only non-zero if ^ ^ 
^ y because the initial cut-off function is constant otherwise. In this 
region x, y) extends smoothly to t — j- and satisfies hl{t, x, y) = 0{t°°) 
by lemma E. Ill 



14 



Finally, if x lies in the Vi^ neighbourhood of a vertex v this property 

3 

follows from lemma 14.111 If x lies in Vi^ on an edge between two vertices v 

3 

and v' this again follows from 14.11] and Xv{x) + Xv'{x) = 1 because this is a 
partition of unity. □ 

5 Construction of the heat kernel 

Definition 5.1. 

The convolution of two continuous kernels P,Q E C*^(]R>o, G, G) on a metric 
graph G is defined as follows. 

{P * Q){t,x,y) := / / P{s,x,z)Q{t ~ s,z,y)dzds 
Jo Jg 

Lemma 5.2. 

Suppose P(t, X, y) = 0{t^) and Q{t, x, y) = 0{t^ ), with k, k' > —1, then 

{P*Q){t,x,y) = 0{t^+^'+') 
Proof. This follows from a direct computation. 

\{P*Q){t,x,y)\^ [ [ Cps^Cgit-sfdzds 
Jo Jg 

k+k'+l / (J\k(^ J\k'^J 



--CpCQCt''+^^' {s'Y{i-s'Ydg 

Jo 

--CpCQCt''+^'+^B{k + l,k' + 1) 



Here B{k + 1, A;' + 1) is the Beta-function and C is the total edge length of 
the graph G. □ 



Definition 5.3. 

Let 



9k{t, X, y) := {dt + Dy)hk{t, x, y) 



Lemma 5.4. 

Let P he a continuous kernel, then 



{dt + Dy){P*hu) = P + P*~gk 
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Proof. Proof adapted from |Ros97j . 

{dt + Dy){P*h)it,x,y) 

= dt / P{s,x,z)hk{t- s,z,y)dzds + Dy{P*hk){t,x,y) 
Jo Jg 

= lim / P{s, z)hk{t — s, z,y)dz + / / P{s,x, z)dthk{t — z,y)dzds 
"^^Jg Jo Jg 

+ / P{s, X, z)Dyhk{t — z,y)dzds 
Jo Jg 

= P{t,x,y) + {P*gk){t,x,y) 
where we used the fact that hk{0, x, y) = Sx{y) in the last step. □ 

Lemma 5.5. 

We have 

oo 

1=1 

in particular, this infinite sum converges. 

Proof. Proof adapted from |Ros97j . By lemma 14.191 

\h{t,x,y)\^Ct'-'2 ^CT'^--^ =:C 
for all ^ t ^ T for some T. We will now show by induction 



rk{t,x,y)\ ^ 



(A:-i + l)(A;-i + 2)...(A;-i + /-l) 
we have just shown the case 1 = 1. We have 
\9*kit,x,y)\ 

^ \9f~^\s,x,z)\-\gk{t-s,z,y)\dzds 



Jo JG 

^ Jo ya(^-| + l)(^-| + 2)...(fc-i + ^-2) 

f-TF^^^^ nl—l ft 



{k-\ + l){k - i + 2) . . . (A; - 1 + / - 2) io 

which finishes the induction. Applying the ratio test shows that the sum of 
the upper bounds converges, this implies that the original sum converges as 
well. □ 
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Theorem 5.6. 

Let 

oo 

e(t, X, y) := hk{t, x, y) + hk{t, x, y) * ^(-l)'^fc(f, x, y)*' 

then e{t, x, y) G C°°(M>o x G x G), is independent of k for k ^ 1 and is the 
heat kernel on G. 

Proof. By lemma 15.51 we know that the infinite sum converges. 
Using lemma 15^ we get: 

oo oo 

{dt + Dy)e{t,x,y) = gk + Y,i-'^y9k{t,x,yy' + gk*Y.^-iygk{t,x,yy' 

1=1 1=1 

= 

We also have 

e{t,x,y)'ilj{y)dy 

« oo 

hk{t,x,y)i){y)dy + / hk*y^{-lfgkit,x,yy^{y)dy 

IG Jg 

^t^o ^(x) + 
The second term doesn't contribute because we have 

oo 

h{t,x,y)*J2{-iy~gl'{t,x,y) = 0{t'') 



G 



by lemmata 15.51 and 15. 2[ □ 

6 The asymptotics of the heat kernel 

Theorem 6.1. 

The heat kernel is approximated by the parametrix. 

e{t,x,y) = hk{t,x,y) + 0{t^) 
Consequently, it admits an asymptotic expansion of the form 



e{t, X, x)dx ~f^.o+ 

G 



oo „ 

y=J2 «"(^)^^ ■ + i E E E (^) ■ ^" 



with some universal coefficients a„(x) and a^(f ) 
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Proof. This first claim follows from * Yl'i^i{~^y 9k' ~ 0{t^) by lemmata 
1531 and [O 

For the second one, we will split the proof into two parts, depending on 
whether x is in the Vj^ neighbourhood of a vertex or not. 

T 

If x is in a V^j^ neighbourhood on the edge a we can unwind the definition 

of the parametrix all the way to the real line. Wer will denote the w-function 
in the coordinates of the star graph around v hj u"". 



e{t, X, x)dx 

-V 



/ hlii^^a,Xa)dXa + 0{t'') 
^ ] / hf, (t, Xa, Xq,) -\- (7^ hj, (t, Xfy, —Xa)dXa -\- 0(t ) 

r'f 1 

)t dXry 

a^v * (=0 

+ /' -^e-^^y"u\{xo^,-x^)t'dx^ + 0{f) 

^ Jo v4:Tct f-f^ 

«~D ^ ^ 1=0 

1 \ - \ - f 3 



^ ^ ^ ^ / (^Xfyj Xa)dXQit 



ar^v 1=0 



k p'-S-t^i 



Now we will treat the part of the integral over , assume that x is on the 

T 

edge a with end vertices v and v\ then we have 

c(tj Xa, Xa'j 

by the properties of the partition of unity x the fact that hk{t,x,y) = 
0{t°°) for X and y bounded away from each other by lemma IT6l Note that 
the part from would be parametrised in the opposite direction but we 
can simply write everything in the same coordinates, this works because the 
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u-functions are local, see 14.41 This means 
e(t, X, x)dx 

Vf 



V / hr{t,x^,x^)dx^ + o{t^) 

^—•LJ'Ji v47rt 7-- 



Putting the two parts together gives 



1 /• 

: / Ui{x,x)dx ■ t'' 

4:7rt ._ J Q 



1=0 



k p'-S-t-i 



+ i E E <° E / ' e-^' -th^)dx^ ■ + 0{t^) 

i)6V ar^v 1=0 



As each ui{x,x) in a coordinate patch is local there is a well defined global 
function ui{x,x), so this implies that ai{x) = ui{x,x). To compute the 
coefficients at the vertices one needs to plug in the Taylor expansions of 
u^{x, —x) for X close to zero. We have 

J2 a'i{v)t' = ^ E <° E / e-^^<(t^x,, -t-2Xo.)dxJ + 0{t^) 

(2) 

We will now compute some w-functions and then use this equation to find 
the first few coefficients. □ 



6.1 The M-functions 

In order to compute the coefficients in the asymptotics, we will compute 
the first few -u-functions explicitly. First, we will need the values on the 
diagonal and second we need the Taylor expansion of ui{x, —x) for 

X approaching zero, that is a vertex of the graph. 
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Lemma 6.2. 

On the diagonal the first terms are 

Uq{x, x) =1 

Ui{x, x) =U (x) 

, , \J"(x) Vixf 

u^{x,x) =^f^^'H^) + \u"{x)U{x) + ^U\xf + \u{xf 

Proof. The w-functions are local and smooth by our assumptions on the po- 
tential. Thus we can carry out this computation globally and coordinate free. 
Assume that x and y are on the same edge. From the transport equation, 
lemma Hl2| we have 

ry 

ui{x,y) = {y - xy^ U{z)dz 

J X 

We will plug in the Taylor expansion of U at the point x. This and all 
following equations are meant to hold for any finite Taylor expansion up to 
a suitable remainder term. 

uii^iy) = {y - X) 2^/ {z-x)dz — - — 
= E(^--)(7TT)! 

This gives the value of ui{x, x) and 

DyUi{x,y) = -}^{y-xy _ 2)1 / + n ~ ^(^) Z^(^ ~ ^ 



(/-2)!(/ + l) ^ (/ + !)! 
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which we now use to compute U2{x, y), vahd for y close to x. 

rv 

U2{x,y) = — {y — x)'"^ / {z — x)DzUi{x, z)dz 

J X 



--{y-x) 



l>2 



I (/-2)!(/ + l) 



+ {y-x) ^ 2^ ^, {z - xr 2^{z - xY+^^^^dz 



(/ + !)! 



|^(/-2)!/(/ + l)^^ 

TT(l)t^\TT(m)f^\ 

\m+l 



+ ^f^J/ + l)!m!(m + / + 2)^^ > 



This determines M2(3;,a^) and we have 



«-4 



[/W(x)[/(-)(x)(m + /)(m + /-l) _ 



Finally, from the transport equation, lemma 14.21 

M3(x, X) =- '^DyU2{x, X) 

\dyU2{x,x) + \u{x)u2{x,x) 



□ 



Lemma 6.3. 
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The Taylor expansions of Ui{x, —x) and U2{x, —x) for small x are given by 

-x) =U{v) + -U"{v)x^ + ^U^^\v)x^ + 0{x^) 
6 5! 

ulix, -x) =-U"{v) + \u{vf + (^l^U^'Kv) + -U{v)U"{v)^ + 0{x') 

where as before we assume that the odd derivatives ofU vanish at the vertices. 

Proof. Recall that in the coordinates of a star graph corresponds to the 
central vertex and all edges are parametrised away from the vertex. Using 
the Taylor expansion of U at zero gives 

u,{x,y) = -{y-x)-'Y.^x^^'-y'^''^'^'^^^^ 



1^0 

^{1 + 1)1^ y 

i>0 ^ ^ m=0 



(/ + !)! 



As in the proof of lemma 1612] this holds for any finite Taylor expansion up to a 
suitable remainder term. This determines the Taylor expansion of ui{x, —x). 

Dymix, y) = -J2 Yl(^-m){l-m- \)x'y'—^ 

1^2 ^ '' m=0 

We will now use this expression to compute the expansion of U2{x, —x) 
for X close to zero. 

1 r 

U2{x,-x)=—J {y - x)DyUi{x,y)dy 

= - i (2/ - E Yl(^-m){l-m- V)x^'—'dy 

=\u"{Q) + ^U^'\^)x^ + \u{^f + lf/(0)t/"(0)x2 + 0{x') 
o oU lb 

□ 
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6.2 Computing the coefficients 

Using the approximation of the heat kernel by the parametrix and the values 
of the M-functions we can now compute the values of the coefficients in the 
asymptotics. 

Theorem 6.4. 

With our assumptions on the potential a^(f) = whenever n is not an in- 
teger. The first few non-zero coefficients are as follows. 

ao{x) = 1 CLoiv) = 1 

ai{x) = U{x) aliv) = U{v) 

a,{x) = lu'\x) + lu{xr aliv) = \u"{v) + ^Uiv)' 

asix) = ^U^'\x) + ]:U"{x)U{x) + ^U'{xY + ]:U{xf 
oU 12 b 

4{v) = j^U^'Hv) + \u{v)U'\v) + ^Uivf 



Proof. We showed ai(x) = ui{x, x) in the proof of theorem 16 -H so the values 
for ao{x) , . . . , a3{x) follow directly from lemma 16^ 

To compute the boundary coefficients we will use the Taylor expansions 
computed in lemma 16.31 and plug them in equation ([2]) from the proof of 
theorem 16.11 Note that because of our assumptions on the potential, only 
the even terms in the Taylor expansions appear, so all terms at half powers 
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of t vanish. 



t — U, 2 1-^: 2 ' 



Here we used 



/ e"" (ix =y- + 0{t 
Jo 2 

Jo 8 



4 ^ 

.'-ot-i 



□ 



Remark 6.5. One can compare the coefficients computed in theorem 16.41 with 
results for manifolds in |Gil79] on the unit interval and the coefficients match. 
The manifold expansion contains boundary coefficients at half powers of t 
but these vanish in our setting because of the assumptions we made on the 
potential. 



24 



Acknowledgement 



This work was made possible through a grant from the DFG, that allowed 
me to spend a year at Royal HoUoway to work on this project. It is a 
pleasure to thank Jens Bolte for coming up with the topic and multiple 
fruitful discussions throughout my work on this. 



References 

[BE08] Jens Bolte and Sebastian Endres, Trace formulae for quantum 
graphs, Analysis on graphs and its applications, Proc. Sympos. 
Pure Math., vol. 77, Amer. Math. Soc, Providence, RI, 2008, 
pp. 247-259. 

[BE09] , The trace formula for quantum graphs with general self 

adjoint boundary conditions, Ann. Henri Poincare 10 (2009), no. 1, 
189-223. 

[BGM71] Marcel Berger, Paul Gauduchon, and Edmond Mazet, Le spectre 
d'une variete riemannienne. Lecture Notes in Mathematics, Vol. 
194, Springer- Verlag, Berlin, 1971. 

[FKW07] S. A. Fulling, P. Kuchment, and J. H. Wilson, Index theorems for 
quantum graphs, J. Phys. A 40 (2007), no. 47, 14165-14180. 

[Gil79] Peter B. Gilkey, Recursion relations and the asymptotic behavior 
of the eigenvalues of the laplacian, Compositio Mathematica 38 
(1979), no. 2, 201-240. 

[Gre71] Peter Greiner, An asymptotic expansion for the heat equation. 
Arch. Rational Mech. Anal. 41 (1971), 163-218. 

[Gri07] Daniel Grieser, Monotone unitary families, 

http://arxiv.org/abs/0711.2869 (2007). 

[KPS07] Vadim Kostrykin, Jiirgen Potthoff, and Robert Schrader, Heat ker- 
nels on metric graphs and a trace formula, Adventures in math- 
ematical physics, Contemp. Math., vol. 447, Amer. Math. Soc, 
Providence, RI, 2007, pp. 175-198. 

[KS99a] V. Kostrykin and R. Schrader, Kirchhoff's rule for quantum wires, 
J. Phys. A 32 (1999), no. 4, 595-630. 



25 



[KS99b] Tsampikos Kottos and Uzy Smilansky, Periodic orbit theory and 
spectral statistics for quantum graphs, Ann. Physics 274 (1999), 
no. 1, 76-124. 

[KS03] Vadim Kostrykin and Robert Schrader, Quantum wires with mag- 
netic fluxes, Comm. Math. Phys. 237 (2003), no. 1-2, 161-179, 
Dedicated to Rudolf Haag. 

[Kuc04] Peter Kuchment, Quantum graphs. I. Some basic structures. Waves 
Random Media 14 (2004), no. 1, S107-S128, Special section on 
quantum graphs. 

[Min53] S. Minakshisundaram, Eigenf unctions on Riemannian manifolds, 
J. Indian Math. Soc. (N.S.) 17 (1953), 159-165. 

[MP49] S. Minakshisundaram and A. Pleijel, Some properties of the ei- 
genfunctions of the Laplace- operator on Riemannian manifolds, 
Canadian J. Math. 1 (1949), 242-256. 

[MS67] H. P. McKean, Jr. and I. M. Singer, Curvature and the eigenvalues 
of the Laplacian, J. Differential Geometry 1 (1967), no. 1, 43-69. 

[PT87] Jiirgen Poschel and Eugene Trubowitz, Inverse spectral theory. 
Pure and Applied Mathematics, vol. 130, Academic Press Inc., 
Boston, MA, 1987. 

[Ros97] Steven Rosenberg, The Laplacian on a Riemannian manifold, Lon- 
don Mathematical Society Student Texts, vol. 31, Cambridge Uni- 
versity Press, Cambridge, 1997, An introduction to analysis on 
manifolds. 

[Rot84] Jean-Pierre Roth, Le spectre du laplacien sur un graphe, Theorie 
du potentiel (Orsay, 1983), Lecture Notes in Math., vol. 1096, 
Springer, Berlin, 1984, pp. 521-539. 

[RSI 2] Ralf Rueckriemen and Uzy Smilansky, Trace formulae for quantum 
graphs with edge potentials, J. Phys. A: Math. Theor. 45 (2012), 
475205. 



26 



